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ABSTRACT 

We analyze the properties of searches devoted to finding planetary transits by observing simple 
stellar systems, such as globular clusters, open clusters, and the Galactic bulge. We develop the 
analytic tools necessary to predict the number of planets that a survey will detect as a function of 
the parameters of the system (age, extinction, distance, richness, mass function), the observational 
' setup (nights observed, bandpass, exposure time, telescope diameter, detector characteristics), site 

properties (seeing, sky background), and the planet properties (frequency, period, and radius). 
. We find that for typical parameters, the detection probability is maximized for /-band observations. 

fSJ ' At fixed planet period and radius, the signal-to-noise ratio of a planetary transit in the /-band is 

weakly dependent on the mass of the primary for sources with flux above the sky background, and 
falls very sharply for sources below sky. Therefore, for typical targets, the number of detectable planets 
is roughly proportional to the number of stars with transiting planets with fluxes above sky (and not 
necessarily the number of sources with photometric error less a given threshold). Furthermore, for 
rising mass functions, the majority of the planets will be detected around sources with fluxes near 
sky. In order to maximize the number of detections, experiments should therefore be tailored such 
| that sources near sky are above the required detection threshold. Once this requirement is met, the 

number of detected planets is relatively weakly dependent on the detection threshold, diameter of 
' the telescope, exposure time, seeing, age of the system, and planet radius, for typical ranges of these 

\Q \ parameters encountered in current transit searches in stellar systems. The number of detected planets 

is a strongly decreasing function of the distance to the system, implying that the nearest, richest 
clusters may prove to be optimal targets. 

■ Subject headings: techniques: photometric - surveys - planetary systems 

o 

^ ■ 1. INTRODUCTION 

Although radial velocity (RV) searches have provided an enormous amount of information about the ensemble 
Q ' properties of extrasolar planets, the interpretation of these results has been somewhat complicated by the fact that 
\ the planets' properties have been shaped by the poorly-understood process of planetary migration. Short-period 
CZ2 . planets (periods P < 10 days, i.e. "Hot Jupiters") are essential for understanding this phenomenon, since they have all 
' almost certainly reached their current positions via migration, and because they are the easiest to detect via several 
methods, including both radial velocities and transits. Thus, it is possible to rapidly acquire the statistics necessary for 
uncovering diagnostic trends in their ensemble properties, which may provide clues to the physical mechanisms that 
drive migration. Although RV searches have been and will continue to be very successful in detecting these planets, 

\ transit searches are rapidly gaining in importance. 

There are currently over a dozen collaborations searching for planets via transits (see lHornell2003|) . These searches 
have recently start ed to come to fruition, and six close-in extrasolar giant planets have been detected using the transit 
techn ique to date <|Konacki et al.ll2003at iBouchv et aHl2"f)0ll iPont et alJl2u0l iKonacki et aTJl2C)0"4l 1207)1 lAlonso et all 



2004), with many more likely to follow. Notably, transit searches have already uncovered a previously unknown 
population of "Very Hot Jupiters" - massive planets with P < 3 days. Current transit searches can be roughly 
divided into two categories. Shallow surveys observe bright (VS. 14) nearby stars with small aperture, large field-of- 
view dedicated instruments llBakos et al.l | 2004tlKane et alJl2004t iBorucki. et al.ll20Q]l iPepper et al.ll2004t lAlonso et alJ 
12004 iMcCullough et al.l2004UDeeg et alJ2004|) . The goal of these surveys is primarily to find trans iting planets around 
bright stars, which facilitate the extensive follow-up studies that are possible for tran siting planets fCh arbonneau et alJ 
l2002HVidal-Madiar et"aLll2003L 120(73 iCharbonneau et alJl2005t ITJeming et alJl2005l) . On the other hand, deep surveys 
monitor faint (V > 14) stars using larger aperture telescopes with small field-of-view instruments. Typically, these 
searches do not use dedicated facilities, and thus are generally limited to campaigns lasting for a few weeks. In contrast 
to the shallow surveys, deep surveys will find planets around stars that are too faint for all but the most rudimentary 
reconnaissance. However, the primary advantage of these searches is that a large number of stars can be simultaneously 
probed for transiting planets. This allows such surveys to detect relatively rare planets, as well as probe planets in 
very different environments, and so robustly constrain the statistics of close-in planets. Deep searches can be f urther 
subd i yided into two categories, n amely searches around field stars in the Galactic plane iUdal skiet al.l2002albl lcl. I2003L 
20041 iMallen-Ornelas et al]l2003[K and searches toward simple stellar systems. 

1 Ohio State University Department of Astronomy, 4055 McPherson Lab, 140 West 18th Ave., Columbus, OH 43210 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 
Electronic address: peppcr@astronomy.ohio-state.edu, sgaudi@cfa.harvard.edu 



2 



CLUSTER TRANSIT SEARCHES 



Simple stellar systems, such as globular clusters, open clusters, and the Galactic bulge, are excellent laboratories 
for transit surveys, as they provide a relatively uniform sample of ~ 10 3-5 stars of the same age, metallicity, and 
distance. Furthermore, such surveys have several important advantages over field surveys. With minimal auxiliary 
observations, stellar systems provide independent estimates for the mass and radius of the target stars through main- 
sequence fitting to color-magnitude diagrams. An independent estimate for the stellar mass and radius, even with a 
crude transit light curve, can allow one to completely characterize the system parameters (assuming a circular orbit 
and a negligible companion mass). Transit data alone, without knowledge of the properties of the host stars, does not 
allow for breaking of the degeneracy between the stellar and planet radius and orbital semi-major axis. As a result, 
considerable additional expenditu r e of resources is requi r ed to confirm the planetary nature of transit candidates from 
field surveys Ipreizler et alj|2002t iKonacki et alJl2003bt iPont et al.ll2005aHBouchv et alJl2005t iGallardo et alJl2005|) . 
Furthermore, using the results of field transit surveys to place constraints on the ensemble properties of close-in planets 
is hampered by a lack of information about the properties of the population of host stars, as well as str ong biases in 
the observed distributions of planetary parameters relat ive to the underlying intrinsic planet population ijGaudi et alJ 
l2004t|Pont et a l. 2005a; Gaudi 2 0051 iDorsher et al]l2005|) . In contrast, the biases encountered in surveys toward stellar 
systems are considerably less severe, and furthermore are easily quantified because the properties of the host stars are 
known. This allows for accurate calibration of the detection efficiency of a particular survey, and so enables robust 
inferences about the population of pl anets from the detecti on (o r lack thereo f) of individual planetary companions 
(Gillil and et alJlSfiofl IWeldrake et alll20o"fil IMocheiska eFalMom IBurke et alJl2005T) . 

There are a nu mber of proje cts devoted to sear ching for transiting pl anets in stellar systems llGill land et al.1 2000: 
I Burke et al.ll2003 iStreet et al]l2003 iBrymtt et al]l2003 iDrake fc Co"okll2004t Ivon Braun et al.ll2005l: IMocheiska etall 
2005; Wcldrak e et al.l l2005: Hid as et al.ll2005(l . These projects have observed or are observing a number of different 
kinds of systems, with various ages, metallicities, and distances, using a variety of observing parameters, such as 
telescope aperture and observing cadence. Although se veral authors have discussed genera l consideratio ns in designing 
and executing optimal surveys toward stellar systems l)Janeslll996t Ivon Braun et al.ll2005t iGaudil 120001 , these studies 
have been somewhat fractured, and primarily qualitative in nature. To date there has been no rigorous, quantitative, 
and comprehensive determination of how the different characteristics of the target system and observing parameters 
affect the number of transiting planets one would expect to find. To this end, here we develop an analytic model of 
transit surveys toward simple, homogeneous stellar systems. This model is useful for understanding the basic properties 
of such surveys, for predicting the yield of a particular survey, as well as for establishing guidelines that observers can 
use to make optimum choices when observing particular targets. 

We concentrate on the simplest model that incorporates the majority of the important features of transit surveys 
toward stellar systems. We consider simple systems containing main-sequence stars of the same age and metallicity. 
We ignore the effects of weather, systematic errors (except at the most rudimentary level), and variations in seeing and 
background. Although we feel our analysis captures the basic properties of such searches without considering these 
effects, it is straightforward to extend our model to include these and other real-world effects. 

In fJ21we develop the equations and overall formalism that we use to characterize the detection probabilities of certain 
planets in specific systems with a given observational setup. In JHwe describe various analytic approximations we use 
to make sense of our detailed calculations, and we show how the transit detection probabilities depend on stellar mass 
and the characteristics of a particular survey. In 21 we list various physical relations and numerical approximations we 
use to calculate detection probabilities. In SjHlwe describe the dependence of the detection probabilities on the input 
parameters, and we present an application of our results in fJ13 We summarize and conclude in <J3 



2. GENERAL FORMALISM 



2.1. The Number of Detected Transiting Planets 

For a given stellar system, the number of transiting planets with periods between P and P + dP and radii between 
r and r + dr that can be detected around stars with masses between M and M + dM is, 

d 3 N det d 2 p dn 

dMd^dP = N * fp drdP rt ° tiM > P > r) dM- (1) 

Here, Nd e t is the number of detected transiting planets; is the total number of stars in the system; d 2 p/drdP is 
the probability that a planet around a star in the system has a period between P and P + dP and a radius between 
r and r + dr; f p is the fraction of stars in the system with planets; Vtot(M, P,r) is the probability that a planet of 
radius r and orbital period P will be detected around a star of mass M; and dn/dM is the mass function of the stars 
in the system. 

There are a number of assumptions that enter into equation Q): 

• We assume that f p and d 2 p/drdP are independent of M. We normalize d 2 p/drdP to unity over a specific range 
of planetary radii and periods, and normalize dn/dM to unity over a specific range of stellar masses. Therefore, 
N* is the number of (single) stars in the mass range of interest, and f p is the fraction of such stars harboring 
planets in the range of planetary radii and periods of interest. The number of such planets is thus N p = f p N*, 
and the fraction that are detected is fdet = Nd e t/N p . The normalization of dn/dM is described in M4.2I and the 
normalization of d 2 p/drdP is described in 
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• We choose to use P as our independent parameter rather than semi-major axis a, since it is the more directly 
observable quantity in transit searches, and it simplifies the following discussion considerably. 

• We note that one of the primary simplifying assumptions in equation Q is that all the target stars are at the 
same distance from the observer, which is an excellent assumption for most stellar systems. 



2.2. Detection Probabilities Vt, "Pw> Vs/n 

Following [Haud3 {2000), we separate VtotiM, P, r) into three factors 

V tot {M, P, r) = VtiM, P)V s/N iM, P, r)7V(P), (2) 

where Vt is the probability that a planet transits its parent star, Vs/n is the probability that, should a transit occur 
during a night of observing, it will yield a Signal-to-Noise ratio (S/N) that is higher than some threshold value, and Vw 
is the window function which describes the probability that more than one transit will occur during the observations. 

2.2.1. Transit Probability Vt 
The probability that a planet will transit its parent star is simply 

1/3 

Vt= *L=( ±L.\ M-^RP- 2 / 3 . (3) 



a \ G 

This form of Vt assumes that the planet is in a circular orbit. We will make this assumption throughout this paper. 



2.2.2. Window Probability Vw 

The window function VwiP) quantifies the probability that a planet with a given period P will exhibit n different 
transits during the times when observations are made. See lGaudT pOOO'l for a mathematical definition of Vw- We will 
consider observational campaigns from single sites comprising a total of N n contiguous nights of length t ni g} lt . For an 
ex ploration of the e f fects o f alternate observing strategies on Vw, we refer the reader to a comprehensive discussion 
bv Ivon Braun et all l)2005f) . We will assume that no time is lost to weather. Finally, we require only that the center 
of the transit occurs during the night; therefore Vw depends only on n, N n , t n i n ht, and P, and does not depen d on 
the transit duration. Note that our definition differs slightly from the definition bv lMall^n-Ornelas et alJ (|20D8T l . In 
Figure QJ, we show Vw as a function of P for N n = 10,20,40 nights and t n i g ht = 7.2 hours, for the requirement of 
n = 2 transits (which we will require throughout). 



2.2.3. Signal-to-Noise Probability Vs/n 
In this section, we determine Vs/n, the probability that a single transit will exceed a S/N value larger than some 
minimum threshold S/N value.* The signal-to-noise of a single transit is S/N = (Ax 2 ) 1 / 2 , where A\ 2 is the difference 
in x" 2 between a constant flux and a transit fit to the data. For simplicity, we will model all transits as boxcar curves. 
In this case, and under the assumption that only a small fraction of the data points occur during transit, the A% 2 of 
a transit is simply, 

A x l = N tr ^ . (4) 

Here N tr is the number of observations during the transit, 5 is the fractional change in the star's brightness during 
the transit, and er is the fractional error of an individual flux measurement. 

The number of observations N tr during a transit is related to the observing timescales: N tr = t tr /it rea d + texp), 
where tt r is the duration of the transit, t rea d is read time of the detector, and t exp is exposure time.§ We can put tt r 
in terms of the fundamental parameters: 

^ = 2x {^ = x {^mT ■ (5) 

$ By folding an observed light curve about the proper period, it is possible to improve the total S/N over that of a single transit by 
~ ra 1 / 2 , where n is the number of transits occurring when observations are made. We have chosen a more conservative approach of requiring 
a minimum S/N based on a single transit because, for observational campaigns such as those typically considered here, the probability of 
seeing many transits is low, and furthermore detailed and well-sampled individual transit signals are crucial for distinguishing bona fide 
transits from false positives. In Appendix ^] we rederive the results of this section for the alternative detection criterion based on the 
total S/N of folded transit light curves. The difference between these two approaches is relatively minor for the surveys considered here, 
although the total S/N approach favors short period planets more heavily. 

§ This model assumes that all transits are observed from beginning to end. We consider the effects of partial transits in Appendix IbI 
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Fig. 1. — The window probability Vw j which is the probability that a planet with a given period will exhibit n separate transits during 
the times observations are made. Here we have assumed 7.2 observable hours each night, that the run is N n nights long, and we have 
required n = 2 transits to occur during the observing window. The lines show the results for N n = 10 (dotted), 20 (solid), and 40 (dashed). 



Here, x is half the length of the chord that traces the path of the transiting planet across the face of the star. 
Geometrically, x = R^Jl — b 2 where b is the impact parameter of the transit. That is, bR is equal to the distance from 
the equator to the latitude of the transit. For a transit with an inclination of 90 degrees, x = R and b — 0, while for a 



grazing eclipse x is nearly and 6=1. 
and therefore t tr = t eq y/l — b 2 . 



We define t eq as the duration of an equatorial transit (i.e. t eq = t tr (b = 0)), 



We assume a transit will be discovered if and only if Axf r is larger than some threshold value Ax„ irl . We note that 



b 2 . 



Therefore, the probability of achieving sufficient S/N is essentially a step function, such that: 
dV S /N 



db 



e 



b 2 - Ax 



min 



(6) 



where is the step function (Q(x) = for x < 0; Q(x) — 1 for x > 0). Equation |JB} provides us with the probability 
that a transit with impact parameter between b and b + db will yield a sufficient S/N to be detected. This can be 
determined for a given set of intrinsic parameters of the system (M, r, and P) and the observational parameters which 
we will list later. 

We can assume that the impact parameters of transiting systems are distributed uniformly. We will take b as our 
fundamental test of S/N, so that if a transit in a system with a given set of intrinsic parameters achieves sufficient 
S/N to be detected with an equatorial transit 6 = 0, then it will be also detectable with any b up to some inclination 
bmax, beyond which point it will not achieve sufficient S/N. We will integrate equation © from b — to b max , which 
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is the range over which Vs/n = 1: 

nl dV s/N 







AxlVl^V 2 - A x 2 mm db. (7) 



This formulation makes it easy to eliminate the step function, since any argument with b > b max will cause the 
argument of the step function to be less than 0, and so the value of the integrand will equal 0. We can therefore 
integrate Vs/n from to b max , so the left hand side of equation (0) becomes 



Vs/n = I db = b max . (8) 
Jo 

We can take the right hand side of equation 10 and note that the argument of the step function will equal when 
evaluated at b max . Setting Ax^l - 6^ - Ax 2 nin = and solving for b max , we then have 



(.eq 



% = 6m«=^-^] , if A Xmm < A xl q , <!)) 



and Vs/n = otherwise. 

We must now determine the dependence of the various factors in equation @ on the independent parameters M, r, 
and P, as well as the observing parameters. Using equation we can put &Xeq m terms of the independent 
parameters and (8/a) 2 . We must therefore relate (S/a) 2 to the independent parameters. Assuming Poisson statistics, 
<y = \/Ns + Nb/Ns, where Ns is the number of photons recorded from a target star in a given exposure, and Nb is 
the number of background photons. In terms of the observing parameters, Ns — f\t exp ir(D/2) 2 , where f\ is the flux 
of photons with wavelength A from the target star, D is the telescope aperture, and we have assumed a filled aperture. 
Flux is related to luminosity by 

h = ^io-^/- (10) 

where L\ is the star's photon luminosity at wavelength A, A\ is the interstellar extinction at wavelength A, and d is 
the distance to the system. Turning to the background sky photons, we can define 

Nb = S skvA nt exp Tr(D/2) 2 , (11) 

where S s ky,X is the photon surface brightness of the sky in wavelength A and f2 is effective area of the seeing disk. 



Putting all this together, we can write Ax 2 Q in terms of the parameters of the planet, primary, and observational 



setup, 



This form can then be inserted into equation © to find Vs/n- 

Note that we have assumed Poisson statistics, no losses due to the atmosphere, telescope, or instrumentation, and 
no additional background flux other than that due to sky (no blending). In we introduce a systematic floor to 
the photometric error a. However, other than this one concession to reality, our results will represent the results of 
ideal, photon-limited experiments, and are therefore in some sense the best case outcomes. When designing actual 
experiments, such real-world complications need to be considered carefully to ensure that they do not substantially 
alter the conclusions drawn here. 



3. ANALYTIC APPROXIMATIONS - SENSITIVITY AS A FUNCTION OF PRIMARY MASS 

To lowest order, the main-sequence population of a coeval, homogeneous stellar system forms a one-parameter 
system of stars. Therefore, a novel aspect of transit searches in stellar systems is that, once the cluster, planet, 
and observational parameters have been specified, the sensitivity of different stars can be characterized by a single 
parameter, namely the stellar mass. This simple behavior, combined with assumptions about the mass-luminosity 
relation, mass-radius relation, and mass function, allows us derive analytic results for the sensitivity of transit surveys 
as a function of stellar mass. 

Here we consider the sensitivity of a given transit search to planets of a given radius r and period P as a function 
of the primary mass M. Adopting power-law forms for the mass-luminosity and mass-radius relations, we rewrite 
the analytic detection probabilities for Vt and Vs/n that we derived in ^2.21 in terms of M. We note that, due to 
the manner in which we have defined it, Vw depends only on P and the observational parameters, and not on M. 
This simplifies the understanding of the sensitivity considerably, since Vw is the only factor that must be calculated 
numerically. 
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R = R, 



© 



(13) 



3.1. Mass-Luminosity and Mass- Radius Relations 
We adopt generic power-law mass- luminosity and mass-radius relations, 

M \ a _ /M x/3a 

Mo J ' iA - LA < Ue, 

where £a,© is the photon luminosity at a wavelength A for a solar-mass star. The power-law index for the mass- 
luminosity relation is wavelength-dependent, such that the f3\ index accounts for bolometric corrections for particular 
bandpasses. 

We note that neither empirically calibrated nor theoretically predicted mass-radius and mass-luminosity relations are 
strict power laws. However, the power-law relations lead to useful analytic results that aid in the intuitive results of the 
more precise results presented later. Furthermore, for stars near M ~ M Q and optical bandpasses, this approximation 
is reasonably accurate. 

For the most part, we will keep the resulting analytic expressions in terms of the variables a and (3\, rather than 
substitute specific values. However, as will become clear, some interesting properties of these expressions are seen for 
realistic values of these parameters. Therefore, where appropriate, we occasionally insert numerical values for a and 
As we show later, for most targets, the J-band proves to be optimal in terms of maximizing the signal-to-noise 
ratio of detected transits. For the /-band, and 0.3 M Q < M < 2M Q , typical values are a = 1 and 0j = 3.5. 



3.2. Dependence ofVs/N on M 
We first consider Ax\ q and Vsin- Substituting equation (|13fl into equations © and l|12jl. we find after some algebra, 




-(3a-/3 A + l/3) 



2(3a-^+l/3) 



1+C 2 



M 

Mr 



1/2 



(14) 

(15) 



where we absorb all the constants and parameters except for mass into the new constants C± and C2 , which are given 

by, 



d = (1024^) 1 / 3 A X L„ ( 1 



r 



D 



GMr. 



1/3 



10 



0AA x 



Co 



4ird 2 S sky ^tt 
L AjQ 10-°- 4 ^ 



(16) 



(17) 



Note that C2 is simply the ratio of the flux in the seeing disk to the flux of a star of M . 

Inspection of the behavior of equation l|15f) as a function of M reveals that there are two different regimes. In the 
first regime the second term within the square brackets is much smaller than unity and hence negligible. This is the 
regime in which the photon noise is dominated by the source (i.e. the target star). In the opposite regime, where 
that term is much larger than unity, the noise is dominated by the sky background. The transition between these two 
regimes occurs at the mass M skv where the flux from the star is equal to the flux from the sky background, 



M, 



sky 



M. 



sky 



If Ay 2 /Ay 2 



(18) 
, < 1 



The behavior of V$/n as a function of mass depends on the value of Axl q / Ax^in & t M 
at M = M s k v , then the ability to detect planets is limited by the source noise for all the stars in the system. Conversely, 
if Axlq/ Axmin > 1 at M = M s k y , then the ability to detect planets around the faintest stars in the system is limited 
by noise due to the sky background. That is to say, a particular experiment can be characterized by whether the flux 
of the faintest star around which a planet can be detected is brighter or dimmer than the sky. We call these the "source 
limited" and "background limited" regimes, respectively. We shall see the implications of this distinction shortly. An 
experiment is in the background limited regime when A% 2 „/A% 



In the source noise limited regime, we find that 



(3a- 



''I 

Px+l/3)/Px 



2 mm > 1 for M > M sky , which implies, 



< 1. 



Ax; 



A X 



2 

rain 



1 

C~i 



where Vs/n becomes 



V 



S/N 



1-Cf 



M 

Mo 

M 

Mq 



-(3q-/3 a +1/3) 



2(3a-/3 A + l/3)' 



1/2 



(Source Limited), 



(Source Limited). 



(19) 



(20) 



(21) 
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On the other hand, in the background noise limited regime, 

L/3) 

(Background Limited), (22) 



A x 2 eq 1 ( M N -(3-^ A+ i/3) 



Ax^ m C X C 2 \M { 
and 



T^S/n 



. 2 , M 



1 - ^ I 



2(3ct-2/3 A + l/3) 



1 1/2 

, (Background Limited). (23) 



Both of these equations have the same general form. For masses below a certain threshold, M t h, there is no chance 
of detecting a transit. The formula for M t h can be determined separately for the two different noise regimes. In the 
source noise limited regime, we have 

M th , a = ( Cl )-V(3a-/Jx+l/3) Mq; (24) 

while in the background noise limited regime, 

M tKb = ( ClC2 )-l/(3a-2/3 x+ l/3) Mq _ (25) 

Thus Vs/n as a function of stellar mass is approximately a step function, and the placement of the step, Mth will 
depend on whether the faintest star around which a planet is detectable (for which A\l q > Ax min ) is brighter or 
dimmer than the sky. Although the labels "source limited" and "background limited" refer to the faintest star for 
which a planet is detectable, and not to all the stars in the system, we shall see shortly that the integrated detection 
probability will depend primarily on the lowest-mass stars. 

It is highly instructive to insert numerical values for a and (3\ and consider the behavior of Ax 2 eq and Vs/n m 
the source and background limited regimes. Adopting values appropriate to the /-band, (a = 1 and (3i = 3.5), we 
have 3a — (3\ + 1/3 = —1/6, and thus A\e q oc M 1//6 in the source-noise limited regime. Thus, for sources above sky, 
the signal-to-noise is an extremely weak function of mass. On the other hand, for sources below sky, we have that 
3a — 2[3\ + 1/3 = —11/3, and thus Ax\ q oc Af 11 / 3 , an extremely strong function of mass. Taken together, these results 
imply that, if it is possible to detect transiting planets around any stars in the target system, it is possible to detect 
planets with the same radius and period around all stars in the system above sky. For stars fainter than sky, the 
detection rapidly becomes impossible with decreasing mass. These effects are illustrated in H5.ll 

These results have an interesting corollary that informs the experimental design. If the experiment is background 
limited (i.e. Ax 2 eq {M 8 k y ) > Ax^ in ), then the minimum stellar mass around which a planet is detectable is Mth,b/MQ = 

(CiC 2 ) 3 / n , whereas in the source limited regime M thtS /M Q = (C1) 6 . Since the constants C\ and C 2 depend on the 
parameters of the target system, the experimental setup, and the observational parameters, these scaling relations 
generally imply that the yield of experiments in the background limited regime is relatively insensitive to the precise 
values of these parameters, whereas the opposite is true for experiments in the source limited regime. Said very 
crudely: specific experiments are either capable of detecting planets or they are not. Experiments should be tailored 
such that Axe ? /AxLn > 1 at M = M sky , which implies that 2C x cf q ^ +1 /3)//3a < ^ but prov id e d this requirement 
is well-satisfied, changing the observational parameters will have little effect on the number of detected planets. 

3.3. Dependence of Vt on M 
We next consider Vt- Substituting equation l|13|) into equation JjJJ and <|12() . 

M \ "~ 1/3 



n = C 3 [ — ) (26) 



where we have defined 

C.3 = „„ ° = 0.238 



\day 



3.4. Dependence of dn/dM on M 

We assume a differential mass function of the form 

dn . / M N 7 



dM^UeJ ' (28) 
The constant k must be chosen such that that the integral over dn/dM is equal to unity, i.e. such that 

dMkl — ) =1, (29) 

M min \ M Qj 

where M max and M m i n are the masses of the largest and smallest stars in the system to be considered. Solving 

f 7 l\M~< +1 - M 7+1 l 

0/ L max mini 



equation ® for k gives us k = (7 + l)Mg/[M^ - Mj +] 
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3.5. Dependence of Vtot{dn/dM) on M 

We can now use these forms for Vs/n an d "Pt, together with assumptions about the mass function of the stellar 
system dn/dM, to evaluate the detection sensitivity to planets with a given set of properties. 

To a first approximation, Vs/n is simply a step function such that Vs/n = Q{M — Mth), where M t h is the minimum 
threshold mass. This is given by Mth = M t h,b if Ax 2 {M s k y ) > A^ in , and Mth = Mth,b otherwise. Thus for masses 
M > Mth, the sensitivity as a function of mass is dominated by the effects of Vt and dn/dM. We can write 

V tot (M,P,r)-^- =Vw(P)C 3 (P)k (-^-) C Q(M-Mt h ). (30) 



dM v ' v ' VM C , 

where ( = a — 1/3 + 7. For a Saltpeter slope of 7 = —2.35, and a = 1, C — —1-68. Therefore, under the assumption 
that the frequency of planets of a given radius and period is independent of the mass of the primary, the number of 
detected planets is dominated by parent stars with mass near M t h, which, in the usual case of a background-dominated 
experiment, is for stars with flux just below the sky. 



4. ADDITIONAL INGREDIENTS 

In fj^jwe adopted several simplifying assumptions and approximations that allowed us to derive analytic expressions 
for the detectability of planets as a function of primary mass. Inspection of these expressions allowed us to infer some 
generic properties of transit searches in stellar systems. However, in order to make realistic estimates of the number 
of planets a particular survey will detect, here we add a few additional ingredients to the basic formalism presented in 
[yj We will also present a somewhat more sophisticated treatment of the mass-luminosity relation, as well as adopt 
specific values for several parameters as necessary to make quantitative predictions. 



4.1. Reconsidering the Mass-Luminosity Relation 

The above analysis approximated the mass-luminosity relation as a simple power law in each wavelength band. As we 
have already discussed, this assumption is incorrect in detail. We therefore provide a somewhat better approximation 
to the mass-luminosity relation. We analytically relate L\ to L& z, assuming purely blackbody emission, and that the 
bolometric mass-luminosity relation can be expressed as a power law: 

M ^ 



Lboi - Lboi.Q \ j ' ( 31 ) 

in which j3 is a single number - the bolometric power law index - instead of the wavelength-depende nt index 0\ i n 
equation l)13|l . Empirically, this is known to be a reasonable approximation for 0.3Mq < M < 2Mq ( Popper 1980). 
We combine this bolometric relation with the mass-radius relation from equation 113(1 . and with Lboi — A-KR 2 aT A . We 
can then write temperature as a function of mass, 

/ M \ (fi-2a)/4 , T \ 1/4 

r < M >= r °(£) • r °=(^i) • (32) 

where Tq = 5777 K is the effective temperature of the sun. We can write the luminosity of a blackbody in a particular 
band X as: 

/+00 
T x {\')B x ,(T){^R 2 ){n)d\' , (33) 
-00 

where B\(T) is the Planck law per unit wavelength, and 7x(\) is the transmission for filter X. We can approximate 
this formula by assuming that the transmission Tx(\) of filter X is a simple top hat with unit height, effective width 
A Ax, and effective wavelength A Cj jc- We can also replace the integral with a product, since B\{T) does not change 
significantly over the intervals defined by the visible or near-infrared filters we will be considering. Also, we can use 
equation (|32|l to write Lx as a function only of mass. Thus, we can rewrite equation (|33|l as 

87r 2 ci? 2 A7t.AA x 

Lx (M) = ^ - . (34) 

eX P \ X C:X kT{M) ) ~ 1 

To check this form of the lum inosity functio n, we compare its reported lumino sities to those from the Yale-Yonsei 
(Y 2 ) isochrones l|Yi et al.ll2001[) . which use the lLeieune. Cuisinier. fc Buserl (|1998j) color calibration. We find that this 
form for Lx(M) is sufficiently accurate for our purposes. In particular, it is much more accurate than the simple 
power-law approximations we considered in <J3 Nevertheless, we find that the qualitative conclusions outlined in that 
section still holds using the more accurate form for the mass-luminosity relation, and thus we can still use the intuition 
gained by studying the behavior predicted by the analytic approximations derived in §^to guide our interpretation of 
the results presented in the rest of the paper. 
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4.2. Normalizing the Mass Function 

To normalize the mass function, we need to determine which values to use for M m i n and M max , which are used to 
compute the normalization constant k. We should choose values that limit the set of stars in the analysis to those 
around which planets are likely to be detected. 

Somewhat anticipating the results from the following sections, we will set M m i n = 0.3Mq for our fiducial calculations. 
This represents the lowest mass star around which a planet can be detected, for typical ranges of the observational, 
system, and planet parameters encountered in current transit searches. In some cases, it may be possible to detect 
planets around stars of lower mass. On the other hand, one might be interested in only those stars for which precise 
radial velocity follow-up is feasible for 8m-class telescopes. Therefore, in M5.4I we consider the effects of varying M m i n 
on the number of detectable planets. 

We set M max to be the most massive main sequence star in the system, i.e. a turnoff star. We determine the mass 
of a turnoff star, M to , using the simple relation Lboi,to = tM to c 2 / A, where e is the net efficiency of hydrogen burning 
(e = 0.00067) and A is the age of the target system. Combining this expression with the bolometric mass-luminosity 
relation from equation l|31ll gives us 

Mto = -j-^-r ■ (35) 

4.3. Minimum Observational Error 

In M2.2.3I we calculate a using a formula for pure photon noise errors. In real observations, photometric errors do not 
get arbitrarily precise for a given source and background. Therefore, we impose a minimum systematic observational 
error of a sys = 0.1% to mimic the practical difficulties of obtaining precise observations of bright stars. The calculated 

/ \ 1/2 

errors therefore become equal to a = ( o~ 2 hot + a 2 ys J , where a p hot is the photon-noise error. 

4.4. Effective Area of the Seeing Disk 

We assume the point-spread function (PSF) is a Gaussian with a full- width half-maximum of 6 see , which has an 
effective area of, 

n =^l d2 see- (36) 

4.5. Saturation Mass 

Detectors have a finite dynamic range, and we clearly cannot detect planets around saturated stars. When integrating 
over mass, we therefore ignore stars with M > M sat , where M sat is the mass of a star that just saturates the detector. 
We assume that a star saturates the detector when the number of photons N p hot from the star and sky that fall into 
the central pixel of the stellar PSF exceeds the full well depth of a pixel, Npw. We approximate N p h ot as, 

- In 2 ' " 



Np ho t = [fx \ 1 - exp 




S sky ,\0 plx t exp TT — , (37) 



where 9 P i X is the angular size of a single pixel. This form assumes a Gaussian PSF perfectly centered on the central 
pixel. The assumption of Gaussian PSF is reasonable for our purposes, and the assumption that the PSF is centered on 
a pixel conservatively underestimates M sat . Formally, equation l|37|) only holds for circular pixels, but is nevertheless 
accurate to < 20% for square pixels. This is sufficient for our purposes. 

4.6. Planet Distribution 

To compute the number of detected planets, we integrate d 3 Nd e t/drdPdM (see equation (QJ) over M, P, and r to 
find Ndet- We therefore must assume a form for the distribution of planets, d 2 p/drdP. We will assume tha t the periods 
are distributed evenly in log space, as is suggested by several analyses fe.g.. lTabachnik fc Tremaine]l2002|) . Since radii 
have been measured for only seven planets, the distribution of the radii is very poorly known. We will therefore simply 
assume a delta function at r = r' , and adopt r' = 0.1i?Q for our fiducial calculations. However, we will also explore 
the detectability as a function of r. Our adopted distribution of periods and radii can therefore be expressed as, 

d 2 n 1 

Hp'a^p^-^' < 38 > 

where AlnP is the logarithmic range of periods of interest. 

From a comparison of the results from radial velocity and transit surveys, it appears that there are two distinct 
populations of close-in massive planets. "Very Hot Jupiters" have periods betw een 1 — 3 days, and are approximately 
ten times less common than "Hot Jupiters" with periods between 3 — 9 days ijGaudi et alJl20 04;'). We will therefore 
consider these two ranges of periods separately. 
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TABLE 1 
Fiducial Parameters 



Parameter Value 



Distance (d) 


2.5 kpc 


Age ( A) 


1 Gyr 


Mass Function Slope (7) 


-2.35 


Bolometric Index 


4.0 


Extinction in I-band (Aj) 


1.25 


Telescope Aperture (D) 


200 cm 


Exposure Time (t exp ) 


60 s 


Seeing (0 see ) 


1 arcsec 


Chi-Square Threshold (Ax* in ) 


30 


Duration of Survey (N n ) 


20 nights 


Planet Radius (r) 


0.1 i?© 


Orbital Period (P) 


2.5 days 



4.7. Extinction 

We consider two models for the extinction. In general, we assume an extinction of a fixed value Aj in the /-band, 
and calculate the extinction in the other bands using the extinction ratios listed in Table [5] We also consider an 
extinction that depends on the distance to the stellar system as, 

Aj{d) = 0.5 mag (j^j , (39) 

where we again use the the extinction ratios in Tabled to determine the extinction in the other bands. We use the 
fixed-extinction law in all calculations and plots unless otherwise specified. 



4.8. Fiducial and Fixed Parameters 

There are a number of parameters in these equations for which we must assign values. In M5.3I we will examine the 
dependence of the detection probabilities on a subset of the most interesting of these parameters. These include the 
cluster distance d, age A, mass function slope 7, and extinction Ai, as well as the telescope aperture D, the exposure 
time t exp , the seeing 6 see , duration of the survey N n , detection threshold Ax„j„, and planet radius r and orbital period 
P. Our choices for the fiducial values of these parameters are listed in Table[Q We do not vary the values of the other 
parameters, either because they are quantities that are empirically well-determined, or because their values are specific 
to the kinds of surveys we are considering here. These quantities are the detector readout time t rea d = 30 seconds, the 
fullwcll depth of the detector Nfw — 10 5 photons, and the angular size of the detector pixels 9 P i X = 0.2 arcsec. We 
assume an exponent of the mass-radius relation of a = 1. We also assume that observations can take place during 7.2 
hours each night, and we require two transits to be observed for a detection. The fiducial values chosen in Tableware 
not intended to represent a specific cluster, but rather to be typical values for star clusters in the Galaxy. 

Since Vs/n depends on the observational band, we calculate Vs/n (and by extension, the overall probability Vtot) 
using 4 different bands, /, V, B, and K, using the \ c x and AAx for each band as defined in lBessell et alJ l|1998|h We 
use the sky brightness in the different bands, /ij, fiy, Hb, Hk, from th e KPNO web s ite.^ Table HI lists the values for 
sky brightness, along with the flux zero point values, which come from lBessell et all l|1998j) . 



5. RESULTS 

We now have all the pieces we need to use equation to evaluate the number of planets Ndet that can be detected by 
a particular survey toward a given stellar system. Our objective in this section is to explore how the overall detection 
probability depends on the various properties of the stellar system, the planets, and the survey, and to provide an 
estimate of the yield of planets for a particular transit survey. 

We begin by exploring the detection sensitivity as a function of host star mass, confirming the basic conclusions 
we derived from our simple analytic considerations presented in SJ3J We then consider the detection probability as 
a function of period, integrated over the mass function of the system. Finally, we consider the fraction of detected 
planets as a function of the various observational and cluster properties, fully integrated over the mass function, as 
well as the assumed planetary period distribution. Unless otherwise stated, we will adopt the fiducial assumptions and 
parameter values described in detail in 21 



" http:/ /www.noao.edu/kpno/manuals/dim/dim.html#ccdtime 
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TABLE 2 

Sky brightnesses, Zero-point Fluxes, 
and Extinction Ratios 



Sky brightness 3 


magnitude per arcsec 2 


W 


20.0 


/"V 


21.8 


MS 


22.7 


MK 


13 


Zero point flux b 


W cm'^nr 1 


/o,2 


1.13 x ID - " 




3.63 x 10" 12 


/o,B 


6.32 x 10" 12 


h,K 


3.96 x 10~ 14 


Extinction ratio c 


Av/Aj 


2.07 


Ab/A! 


2.74 


Ak/Ai 


0.232 



a from the KPNO website, 

http:/ /www.noao.edu/kpno/manuals/dim/dim.html#ccdtimc 



from Bcsscll ct al. ( 1998) 
c from EieY fc Mcrrificld 1: 199811 ■ recalculated 
for this paper using I as the reference band. 



5.1. Sensitivity to Host Star Mass M 

We first consider the sensitivity as a function of the host star mass. We begin by considering Ax^ q versus mass 
for our fiducial parameters and r = 0.1 J?© and P = 2.5 days. This is shown in Figure for the four photometric 
bands we consider. We also show our fiducial value of A\min = 30, and the quantities Mth,i and M s ky,i introduced 
in §'3.2\ In order to elucidate the effects of systematic errors, we show AXe q for our no systematic error, and for our 
fiducial assumption of a systematic error of a sys — 0.1%. When the systematic error is negligible, we find that A\e q 
is approximately independent of mass for M > M s f. y , as anticipated in §'<S.2\ However, when systematic errors are 
included, A\1 q nas a peak, which for our adopted values is near Mq. We also see that, for all of the photometric 
bandpasses and fiducial parameter values we consider, the surveys are in the background-limited regime, and that the 
S/N is highest in the /-band, implying that, all else equal, the number of detected planets will be maximized when 
using this band for these fiducial parameter values. 

It is interesting to note in Figure [21 that the behavior of A\e q versus mass is fundamentally different in K than 
the optical bandpasses. The basic reason for this is that, for the mass range considered here (0.1 < M/M < 2), 
observations in K sample the stellar spectrum in the Rayleigh- Jeans tail, whereas observations in the optical sample 
near the blackbody peak or in the Wein exponential tail. Therefore, the Axlq falls more gradually toward lower masses 
for observations in K. We will see this fundamentally different behavior in K exhibited in many of the following results. 

We next consider the overall detection probability Vtot(M, P, r), and its various components, Vw, Vt, and Vs/n- "Pt 
is described by equations (ffi>|) and l(T7|) : Ps/n is described by equations JTSJ, (QSJ, and lfT7|) : and the window function 
Vw is shown in Figure Q}. We plot these various detection probabilities, and Vtot, versus host star mass for our set 
of fiducial parameter values in Figure 13k- The overall shape of the Vs/n curve in the top plot is simple; it shows that 
a transit will be detected if the star's mass is greater than M t h- The small downturn at high masses is due to the 
systematic error introduced in That is because as M increases, the depth of the transit decreases (because of the 
increasing R) but the photometric precision also increases. However, by placing a limit on the measured precision, at 
high masses the decrease in 5 is no longer offset by a decrease in a, and so the sensitivity dips. 

We combine these pieces with the mass function dn/dM in Figure^. There are a couple interesting features of the 
lower plot. The mass function cuts off a little over 2M Q because that is the turnoff mass for a system with the fiducial 
parameters we are using. The probability curve for K band cuts off before that point, though. That is because a 
detector with the fiducial values we have chosen (D = 2m, t exp = 60s, 8 P i X = 0.2 arcsec, and Npw = 10 5 photons) 
saturates at that mass in K, while the values for M sat for /, V, and B are higher than M to for this fiducial stellar 
system. 

Looking at the lower plot, it is clear that / band is the best one to use to detect planets. The number of stars 
increases with decreasing mass, and it is possible to detect planets around stars of lower mass in the J-band. Since 
Ndet involves the integral of Vtot(dn / dM) over mass, the total number of planets detected will be larger in the /-band. 
We will see in M5. 31 that the / band remains optimal for most parameter combinations encountered in current transit 
surveys. 



5.2. Sensitivity to Period 
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Fig. 2. — ^Xeq versus host star mass for our fiducial parameters and r = O.IRq and P = 2.5 days. The vertical lines show Mth,l 
and M a k y j, while the horizontal dashed line shows our fiducial value for AxJ, in . The dotted lines show the curves of A\e q without the 
inclusion of the systematic error a S ys = 0.1%. 



We next examine the sensitivity as a function of period. We consider the total detection probability Vtot, weighted 
by the mass function, dn/dM, integrated over period, i.e. J r P t ot(dn/dM)dM . This is shown in Figure^] 

The strong sensitivity to shorter period planets is clear, and arises from competition from several effects. The signal- 
to-noise probability Vs/n increases for increasing P, since a planet with a longer period will have a transit with a 
longer duration, and so there will be more observations during the transit and hence higher S/N. However, this effect is 
more than compensated by the fact that the transit probability is oc P~ 2 / 3 , and the window probability Vw generally 
increases for smaller periods (see Figure^), since there is a greater chance of detecting two transits for shorter periods. 

5.3. Sensitivity to Parameters 

In this section, we examine how the fraction of detected planets fdet depends on the various input parameters 
considered and listed in Table Conceptually, there arc three different classes of parameters in Tabled Five of the 
parameters describe the properties of the target system: d, A, 7, /3, and Ai. Five of the parameters are properties 
of the observing setup: D, t exp , 9 see , Ax^ lin , and N n . The two remaining parameters, P and r, are properties of 
individual planets. 
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Fig. 3. — Detection probabilities versus mass. Panel (a) shows Vw, Vt> and Vg/tfy and the product of all three, Vtot- The curves 
for Vtot are displayed with dotted lines, to make them more easily distinguishable from the curves for Vs/n- Panel (b) shows the mass 
function dN/dM, along with the product of it and Vtot- In both plots, Vs/n an d Vtot are shown in four different bands, / (black), V 
(purple), B (blue), and K (red). The plots are calculated using the fiducial parameter values in Table^ and a period of P = 2.5 days, for 
which Vw = 0.63. The cutoff in (b) is due to the fact that for a system with the fiducial parameter values used for this plot, the turnoff 
mass is a little over 2Mq . 



Integrating over mass, period, and radius, the fraction of planets detected is, 

^ = Wrj r 1 7 drdPdM ^ tot (M, P, *) *L . (40) 

In Figure we plot fdet versus planet radius. We plot fdet for the six parameters of the target system in Figure 
El and the five observing parameters in Figure [7] We shall now go though each of the parameters and describe the 
dependencies. 

• Radius - The dependence on planetary radius r shows that detection probabilities increase very quickly up to 
the fiducial value of r = O.li?©, at which point a transit with a planet of that radius has sufficient signal to be 
detected around nearly all the stars in the system. In this plot we also see that while the curves for /, V, and 
B bands all have this similar "step-function" shape, whereas the rise for the if -band is more gradual. 

• Age - From Figure Et, we see fdet is quite insensitive to the age of the system in / band, and is somewhat more 
sensitive in the other bands. This is because, as we see in FigureOt, a larger proportion of the planets detected 
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Period (days) 

Fig. 4. — The total detection probability Vtot, weighted by the mass function dn/dM , and integrated over mass, i.e. / Vtotidnj dM)dM , 
versus planet period, for our fiducial parameters (see Table HI . The window function used for this calculation is the same Vw as m Figure 
111 The various colors represent the different observing bands, I (black), V (purple), B (blue), and K (red). 

in the other bands are at higher masses. The fact that M sat is lower than M to in K accounts for the break in 
the curve in K. 

• Bolometric Index - Figure |SJd shows that fdet is weakly dependent on the value of j3. Therefore our choice of 
(3 — 4.0 is not so important, as fdet is essentially the same for 3.5 < (3 < 5.0, which encompasses the whole range 
of values that are typically used for the mass-bolometric luminosity relation. 

• Distance - This is a key parameter. A nearby system will have many saturated stars, which accounts for the 
turnover at small distances. For fixed extinction (Figure|BJj) , as the system gets further away, the signal drops and 
planets cannot be detected around the smaller stars in the system. For distance-dependent extinction (Figure 
EH), that effect is compounded at large distances, although the extinction has much less of an effect in K. In 
both cases, for sufficiently large distances, the system transitions into the source-limited regime, at which point 
fdet drops precipitously. 

• Extinction - Figure plots the dependence on the value of Aj, showing the effects of greater extinction at a 
fixed distance. In a sense, combining the effects from Figure 0: and Figure gives us Figure 01, although the 
combination is more complex than a simple multiplication. 

• Mass Function - The slope of the mass function 7 determines th e relative numb er of smaller stars and larger 
stars. The fiducial value of 7 — —2.35 is the usual Salpeter slope l)Salpeterlll955j) . As we see in Figure Ef, the 
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Planet Radius (R,) 

Fig. 5. — The fraction of planets detected, f^ et versus planet radius, for our fiducial parameters (see TablcQ. In order to convert to the 
number of planets detected N^ et , these numbers should be multiplied by the total number of stars TV* in the system with masses between 
M = 0.3Mq and the turn-off mass, and by the fraction f p of these stars with planets, i.e. iV<j e j = fdetfpN*. The colors represent the 
different observing bands, I (black), V (purple), B (blue), and K (red). The solid lines show f^ et for "Very Hot Jupiters," i.e. planets 
with periods P = 1 — 3 days. The dotted lines show fdet f° r "Hot Jupiters," planets with periods P = 3 — 9 days. The vertical green line 
indicates the fiducial value we use for the other plots, r = O.IR© ~ IRj. 



detection probabilities do not depend greatly on the exact value of the slope, although for larger values of 7 bluer 
bandpasses become more competitive, as expected. 

• Nights Observed - In Figure [7^, we see that the an observing campaign lasting about 15 nights will detect 
two or more transits from nearly all the "Very Hot Jupiters" (1 day < P < 3 days) that satisfy the detection 
threshold, but to detect two or more transits from most of the detectable "Hot Jupiters" (3 days < P < 9 days), 
the survey should last more than twice as long. Since we assume perfect weather in this analysis, even more time 
should be expected to fully detect the most possible transits. 

• Seeing - An increase in the seeing means an increase in the size of the PSF, and so an increase the number of 
pixels over which the flux of the stars is distributed. This affects fdet is two distinct and opposite ways. First, 
this increases the contribution of the background noise at fixed mass, therefore increasing M t h- Second, this 
decreases the number of photons in the central pixel, and so increases the mass at which the detector saturates, 
M sat . As discussed in H3.21 Mth is rather weakly dependent on seeing for experiments in the background limited 
regime, due primarily to the fact that the mass-luminosity relation is so steep. Furthermore, the increase in M t h 
is partially compensated for by the increase in M sat . As a result, the fdet varies very little for the typical range 
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Fig. 6. — The fraction of planets detected, fdeti as a function of the parameters of the target stellar system. Line types are the same as 
in Figure |5] The vertical green line in each plot indicates the fiducial value for that parameter used in calculating all the other plots. The 
panels show /^ej versus (a) age of the system, (b) power-law index of the mass-bolometric luminosity relation, (c) distance to the system, 
for a fixed extinction, (d) distance to the system, for a distance-dependent extinction, (e) /-band extinction, and (f) the index of the mass 
function. 



of seeing encountered in real observations, as seen in Figure 03. Since the sky is so much brighter in K, the 
seeing dependence is somewhat greater in that band. 

• Exposure Time - There are three effects of t exp . First, a longer exposure time increases the total number of 
photons in a single observation. Second, longer exposure times decrease the number of observations per transit. 
These two effects effectively cancel when t exp 3> tread- The third effect of t exp is that very long exposure times 
cause bright stars to saturate. In Figure 0: the saturation effect is the reason why fdet in K falls so quickly, 
since for our fiducial setup the large number of sky photons alone already brings the pixels close to saturation in 
K . In the other bands complete saturation does not occur even at t exp = 10 minutes. Since saturation involves 
both t exp and D, we shall see in IHJthat the simultaneous consideration of both factors is important. 

• Telescope Aperture - This factor enters in two ways. Larger apertures allow a survey to reach the detection 
threshold for fainter stars, yet also lead to saturation of brighter stars. In Figure 01 in K we see that fdet 
plummets a little past D — 200 cm, since at that aperture the sky photons alone saturate the pixels. In /, the 
situation is complicated. Increasing D decreases M t h- However, looking back at Figure we see that in /, M t h 
is just a little larger than 0.3M Q , which we take as the minimum observable mass. Thus increasing D eventually 
pushes M t h below 0.3M Q . Any further increase in D will therefore not result in additional detections at the 
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Fig. 7. — Same as Fig. El but with f det as a function of the observational parameters. The panels show / £ j et versus (a) number of nights of 
the observational campaign, (b) full-width half- maximum of the point-spread function (the seeing), (c) the exposure time, (d) the diameter 
of the telescope, and (e) the minimum Xmi n required for detection. 



low-mass end, and instead simply lowers fdet as an increasing number of high-mass stars saturate the detector. 
In V and B, Mth continues to decrease for larger apertures, but due to its weak dependence on D, it is always 
above O.3M for D < 2000 cm. Further, the increase in fdet due to decreasing M t h is compensated for by the 
decrease in M satl such that fdet is nearly independent of the aperture for D > 500 cm and B— and V— band 
observations. 

• Detection threshold - The choice of Ax^ in is strongly related to how much follow-up time and resources are 
available for confirming transit candidates. Since false positives are a big hurdle in confirming transits, it is best 
to choose a high value for Ax^ lin - As anticipated in H3.2I and seen in Figure [7^, the dependence of fdet on 
Axf nin is relatively weak until the background-limited regime is reached, at which point fdet falls rapidly. For 
our fiducial parameter values which arc representative of many open cluster surveys, rather stringent detection 
criteria of Ax 2 min ^ 100 can be tolerated without an unacceptably large reduction in the detection efficiency. 



5.4. Minimum Mass 

Up to this point, our results have been normalized such that fdet is the fraction of the planets orbiting stars with 
masses between 0.3 Mq < M < M to that are detected. We have not considered masses below 0.3 M©. The lower 
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mass limit was chosen because this is approximately the minimum mass around which a planet was detectable in the 
/-band for our fiducial assumptions (see Fig. [2J ■ Furthermore, it is also approximately the completeness limit of the 
deepest mass function determinations for rich old open clusters (e.g.. iKalirai et al"]l2001|h 

In some instances, it may be possible to detect planets around stars with masses considerably smaller than we have 
considered, with M < 0.3M Q . Since constraints on planets orbiting such very low mass stars are meager, we briefly 
consider the detectability of planets around host star masses in this regime. Specifically, we perform the same analysis, 
except now we consider stars with masses in the full range Mhb < M < M to , where Mhb = 0.08 M is the mass 
at the hydrogen-burning limit. In order to make these results directly comparable to our previous results, we will 
continue to normalize the mass function such that iV* is the total number of stars between 0.3 M@ < M < M to , and 
Ndet = fpfdetN*, with f p the fraction of stars with planets, and fdet the number of planets orbiting stars between 
0.3 M < M < M to that are detected. In this way, fdet can now formally exceed unity, although in practice this is never 
the case. Figure |Ht» shows fdet versus distance including stars down to the hydrogen-burning limit. We see that, for 
monotonically rising mass functions, it may be possible to increase the number of detections significantly by considering 
very low mass primaries. However, initial mass functions are observed to have breaks near M ~ (0.3 — 0.5) M@, such 
that this boost is probably not realized in practice, and furthermore any detections around such low-mass primaries will 
be quite difficult to confirm, as we discuss below. Nevertheless, the potential for constraining the planetary population 
of very low-mass primaries is noteworthy. 

In order to deter mine planet masses, as well as eliminate the many kinds of a strophysical false positives that mimic 
planetary transits {Torres et aTll2004t iMandushev et"aT1l2005t IPont et alJl2005alK reasonably precise ~50ms 1 radial 
velocity (RV) follow-up measurements of candidate transits are required. Since the majority of the stars probed by 
transit surveys toward stellar systems are relatively faint, the ability to perform RV follow-up to this precision is a 
serious concern. The current state-of-the-art RV measurements on faint stars usin g lOm -class teles copes c an reach 
precisions of - 50 m s" 1 on stars with V < 17 (|Konacki et al.ll2003bUBouchv et al.l l2005: Po nt et al.ll2005H) . It may 
be possible to push this limit to somewhat fainter stars with more ambitious allocation of resources, or improvements 
in future technology. 

In order to estimate what fraction of detected planets can be confirmed using RV follow-up, in Figure|Sl we plot fdet 
versus distance, where we consider only those host stars with apparent U-magnitudes of V < 17 and V < 18. These 
can be compared directly to the case where we consider all stars with M > 0.3 Mq. The first feature of Figure 03 that 
is noticed is that the plots of fdet drop faster for the magnitude limited cases than for the mass limited case, because 
of the nature of the mass-luminosity relation. For our fiducial parameters, it is clear that the number of candidates 
that can be confirmed is considerably smaller than the number that can be detected. Furthermore, the advantage of 
observations in the /-band is effectively removed, since most of the additional /-band detections are too faint for RV 
confirmation. 

It is clear that the ability to perform RV follow-up on candidate planetary transits must be considered carefully 
when designing a transit survey. The question of how to devise a photometric survey that maximizes the number 
of detected planets while accounting for the ability to perform spectroscopic follow-up is outside of the scope of this 
paper, but the formalism we have introduced here should provide the tools to do so. 

6. AN APPLICATION 

The most obvious application of our results is to use the predictions for the number of detectable planets to choose 
optimal targets for a particular survey, and to derive strategies to optimize the number of detected transits for a 
specific target. Since the specifics of the optimal strategies will depend on detailed properties of the survey, such as the 
site, detector, telescope, and time allotment, here we will not attempt a comprehensive discussion, but rather simply 
suggest heuristic guidelines motivated by a couple of specific examples. 

One important question is which targets are optimal in the sense of allowing the largest number of possible detections. 
There are a number of factors that may enter into target selection, including visibility, metallicity, richness, extinction, 
size and distance. We illustrate how our results can be used to quantify and optimize cluster selection, using the 
example of the trade-off between cluster distance and exposure time. As we have shown, fdet is a strong decreasing 
function of cluster distance, such that closer clusters are generally preferred. However, saturation of bright stars is also 
more problematic for more nearby clusters. This can be partially compensated for by decreasing the exposure time, 
but only until t exp becomes comparable to t rea d- In Figure Eh we show how fdet varies with distance to the target 
for various exposure times. From this result, we see that clusters with distances d ~ 2 kpc are optimal, for telescope 
apertures of D ~ 200 cm. 

A somewhat different problem is to determine, given a particular target system in which one wants to search for 
planets, what is the optimal observational setup. If the intent of a survey is to detect all the transits of the brightest 
stars, then the exposure time should be set such that the survey saturates at the turnoff stars of the target system. 
In that case, the exposure time can be calculated using equations (|10|) . (|13|l . 1)37(1 . and H35|l. On the other hand, if the 
intent is to configure the parameters to achieve the largest number of photometric detections, there are two factors 
at which we should look more closely: aperture size and exposure time. We can see in Figure EJi how fdet varies 
with aperture size for various exposure times. When the aperture size is large, saturation effects reduce the detection 
efficiency very quickly, and this can only be partially compensated for by decreasing the exposure time. As a result, 
for transit surveys aiming to detect Jupiter-sized planets, telescopes with apertures of D ~ (200 — 400) cm are optimal. 
Exposure times of less then a couple minutes are generally sufficient. 
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Fig. 8. — Plots of f^, t versus distance. The different plots use different criteria for the minimum mass cutoff, with a fixed mass of 0.3Mq 
in Panel (a), O.O8M0 in Panel (b), or minimum masses corresponding to magnitude cuts V < 17 in Panel (c) and V < 18 in Panel (d). 



These preliminary calculations provide some guide to the best places and optimal methods to look for planetary 
transits. Observers searching for transits can use the formalism derived here to precisely determine which systems to 
search for transits, and what observing setup to use. 

7. SUMMARY AND DISCUSSION 

In this paper we have developed a formalism to predict the efficiency of searches for transits in stellar systems. We 
have taken into account most relevant parameters that affect the number of transits that can be observed, and we have 
described how the total number of expected transit detections depends on these parameters. Our primary results are 
as follows. 

(1) /-band is optimal. For the range of parameters encountered in most transit surveys, observations in the I- 
band maximize the number of detected planetary transits. In general, redder bands are preferred also because 
the effects of limb-darkening are minimized, which aids in the interpretation of transit candidates and in the 
elimination of false positives ( Mallc mOrnelas et al]l2003f) . However we have not taken into account the variation 
in quantum efficiency of detectors as a function of wavelength. This is likely a small effect. For some detectors, 
fringing can be a serious problem in the /-band. Thus, in some cases, somewhat bluer bands (i.e. the i?-band) 
may be preferable. Surprisingly, in essentially no case do we find that observations in the /C-band outperform 
those in the /-band. 
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Fig. 9. — The fraction of detected planets fd et as a function of telescope aperture and distance for various exposure times, showing 
only the I-band curves. The solid lines represent Very Hot Jupiters (1 day < P < 3 days), while the dotted lines represent Hot Jupiters 
(3 days < P < 9 days). 



(2) The S/N depends weakly on primary mass. For the /-band, and assuming only Poisson uncertainties, the 
signal-to-noise ratio of a planetary transit is an extremely weak function of mass for sources with flux above that 
of the sky background. For sources with flux below sky, the S/N is a very strong decreasing function of mass. 

(3) The number of detections is proportional to the number of stars above sky. As a direct consequence 
of item (2), if one can find planets around any stars in the target system, one can detect planets around all 
stars in the system with fluxes above sky. Therefore, the number of planets that are detectable is proportional 
to the number of stars with flux above sky. This is quite distinct from the usual assumption that the number 
of detectable planets is proportional to the number of stars with photometric error less than a given precision, 
usually taken to be ~ 1%. Estimates based on this canonical criteria will typically be incorrect. 

(4) Most planets will be detected around stars with flux near sky. Under the typically valid assumption of 
a mass function that rises toward lower-mass stars, item (3) implies that most planets will be detected around 
stars that have fluxes approximately equal to the flux of the sky background. 

(5) Planets orbiting stars near sky must be detectable. The primary requirement for a successful transit survey 

is that the planets orbiting stars with fluxes near the sky background must be detectable. This requirement is 
formulated mathematically in equation (|19|l . Provided this requirement is met, the number of detected planets 
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is a rather weak function of the radius of the planet, index of the bolometric mass-luminosity relation, age of the 
system, index of the mass function, seeing, exposure time, telescope aperture, and the detection threshold. 

(6) The richest, closest systems are optimal. The number of detected planets has the strongest dependence on 
the distance modulus, the distance and extinction to the cluster. Systems at distances of d < 2 kpc are optimal. 
Very nearby (d < 1 kpc) systems may have difficulties with saturation of bright stars, as well fitting within the 
field-of-view of the detector. 

(7) Follow-up of the majority of candidates may be difficult. The majority of planets in typical target 
systems are likely to be detected around stars with apparent magnitudes of V > 17, making precision RV follow- 
up difficult. This is an important conclusion which may affect the design of transit surveys. Difficulties with RV 
follow-up are partially ameliorated by the fact that surveys toward stellar systems are much less prone to the 
ambiguities with the interpretation of candidate detections encountered in field surveys, since the properties of 
the primaries are better known. 

There are a number of ways in which our analysis could be expanded and refined. For instance, we do not take into 
account the met allicity of the observed s ystems. Studies have indicated that planets are more common around high- 
metallicity stars l)Fischer fc Vale nti 2005) , and as that correlation becomes better characterized we can add metallicity 
to the parameters we examine. It would also be useful to include the effects of observability on the window function, 
which is important for the selection of optimal targets. Our analysis also does not account for bad weather. Further 
work could examine what kinds of inclement weather are most damaging for a transit search and could possibly address 
the question of whether it is possible to partially compensate for inclement weather by adopting more sophisticated 
observing strategies. Another potential refinement would be to account for different forms for the mass function, rather 
than relying on the simple Salpeter shape we use in this paper, such as a broken power-law function with different 
power-law indices for high mass and low mass stars. Lastly, we do not account for stellar binarity, which also generally 
decreases detection probability. 



We would like to thank Chris Burke for useful discussions. We would also like to thank the referee for a prompt 
response and helpful suggestions. This work was supported by a Menzel Fellowship from the Harvard College Ob- 
servatory, and also by the National Aeronautics and Space Administration under Grant No. NNG04GO70G issued 
through the Origins of Solar Systems program. 

APPENDIX 

TOTAL S/N FORMULATION 

In the main text, we derived expressions for the signal-to-noise ratio of a transiting planet S/N = (Ax 2 ) 1 / 2 for a 
single transit. The probability Vs/n that a planet will have a S/N that exceeds a given threshold, and all subsequent 
calculations, were based on this single-transit S/N criterion. However, planets will generally exhibit multiple transits, 
and it is possible, by folding an observed light curve about the proper period, to improve the total S/N over that 
of a single transit by ~ n 1 / 2 , where n is the number of observed transits. In fact, popular trans it search algorithms 
operate on phas e-fold ed light curves, and so trigger based on this total signal-to-noise ratio ijKovacs et alJ 1200*21: 
lAigrain fc Ir^lsffiol IWeldrake fc SackettJ l2005) . It is therefore interesting to rederive our expressions based on this 
total S/N formulation. 

The general expression for the number of detected planets N^et remains the same, but the expression for the total 
detection probability Pt t(M, P, r) needs to be altered, 

V tot (M, P, r) = V t (M, PyPff N (M, P, r)V w (P), (Al) 

where Vt and Vw are the transit and window probabilities as before, and V l gJ N is now the probability that the total 
signal-to-noise ratio is higher than some threshold value. 
The total signal-to-noise probability, "P^v can ^ e derived in an analogous way as the one-transit signal-to-noise 

probability (see H2.2.3|l . We begin by defining dV s ° l N /db = O [Ax 2 r — Axf nin ] , where Ax 2 r is the difference in x 2 
between a constant flux and transit fit to the data, 

A X 1 = N^(^j\ (A2) 

Here iV'°* is the total number of observations taken during any transit, and <5 and a are as before. 

For no aliasing, and periods much shorter than the length of the observational campaign, the total number of 
observations during transit is simply the transit duty cycle itr/P, times the total number of observations Nt t, 



(A3) 
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In fact, for campaigns of finite durations from single sites, aliasing cannot be ignored, and there will be a dispersion in 
the fraction of points during transit about the naive estimate ttr/P- For long campa igns lasting more t han ~ 40 days, 
aliasing effects are generally not dominant, although they are still significant fsee lGaudi et al.l 120041 for examples). 
They can be accounted for by integrating d 2 V '5 i ^ / 'dipdb over the transit phase <fi as well as impact parameter b. For 
simplicity, we will ignore aliasing effects here, and assume equation i|A3|l . 
Since ttr = Vl — b 2 t eq , we can write, 

A x l = A x l q VT^¥, A X l q = N t Jf(t) 2 ,. (A4) 



The total signal-to-noise probability is then just the integral over impact parameter, 

ytot f 1 d ? Sl 

s/N -.L db 

which yields 



n% - t d ^b, ( A5) 



2 



'AX 2 

>S/N = y-\ 

if A X 2 min < A X 2 eq , and T S /N = otherwise. 

We can write \\ q m more explicit terms using the expressions for er, 5, and t eq derived previously. The additional 
new ingredient is the expression for Ntot- If we assume that the campaign lasts N n nights, each with a duration of 
tnight, an d that observations are made continuously, then the total number of data points is 

N tot = K l ht N n . (A7) 

Combining this with the expressions we derived in ^2.2.31 we arrive at the expression, 

which can be com pared to the analogous expression for a single transit, equation 112(1. Comparison of equation l|12|l 
and equation (|A8(I reveals that the ratio of x\ q f° r the total signal-to- noise ratio formulation to \\ q f° r the single-transit 
formulation is cx P _1 . Thus the total S/N formulation favors short-period planets more heavily than the single-transit 
formulation. 



EFFECT OF PARTIAL TRANSITS 



Let us return for the moment to our definition of N tr . This variable represents the number of observations of the 
system during a single transit. We stated earlier that N tr — t tr /(t reae i + t exp ). However, that formula is only valid if 
the entire transit is observed during the night; it does not hold if only partial transits are observed, i.e. if the transit 
begins before the start of the night or ends after the end of the night. In those cases the transit is observed for a time 
less than t tr , the number of observations during transit is less than N tr , and therefore the signal-to-noise is less than 
the naive estimate in ^12.2.31 

To account for partial transits, we rewrite the transit duration as, 



ttr =teqy/l- IP f {</>), (Bl) 

where /(</>) is the fraction of the total transit duration that occurs during the observation window, as a function of 
the phase <fi of the transit. For uniform sampling, and t eq < t n ighu this is simply, 



2 
1 

_L 

2 

I 



tnight 



if 

if 
if 

otherwise 




(B2) 



> 1 



where <f> = is the beginning of the night and <p = 1 is the end of the night. Note also that we have also conservatively 
assumed that a transit cannot be detected if it is observed for less than half of its total duration. 
Following the discussion in H2.2.3I we write 



d 2 V S /N 
dbd(f> 



e 



AxLVl - & 2 /(0) - Ax 
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win 
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0.2 0.4 0.6 0.8 1 

Ay 2 . /Ay 2 

Fig. BIO. — The probability Vs/n that a planet, producing A\ 2 = ^Xeq f° r an equatorial transit, will yield a A\ 2 greater than a given 
threshold AxJ, in , when integrated over all impact parameters and phases, for various values of the ratio of the equatorial transit duration 
t eq to the duration of the night t n i ff ht. The solid line shows the approximation Vs/n = y 1 — (^Xmin / '^xl 9 ) 2 used in the main test, 
which is valid for t eq -C t-night 



Proceeding in the same way as in > 12.2.3I we integrate equation l|B3(l over b from to b m ax, and <f> from to 1, assuming 
a uniform distribution for b and <fi, solve for b max , i.e., 



dbd(j) 



We do not attempt to solve equation l|B4|) analytically, rather we evaluate it numerically, noting that Vs/n depends 
only on the ratios ^Xmin/xlq an( A t eq /t nig ht- Figure lBTOl shows Vs/n as a function of AXmin/xlq f° r equatorial transit 
durations lasting 10— 50% of the night. We also show the result for the simplified assumption of t eq <C tnight that we 
adopted throughout. We conclude that our simple assumption is sufficient for purposes, but note that it overestimates 
Vs/n by as much as 25% for certain combinations of parameters. 
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